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Abstract 

Remodelling of biological tissue, due to changes in microstructure, is treated in the continuum me- 
chanical setting. Microstructural change is expressed as an evolution of the reference configuration. This 
evolution is expressed as a point-to-point map from the reference configuration to a remodelled con- 
figuration. A "preferred" change in configuration is considered in the form of a globally incompatible 
tangent map. This field could be experimentally determined, or specified from other insight. Issues of 
global compatibility and evolution equations for the resulting configurations are addressed. It is hypoth- 
esized that the tissue reaches local equilibrium with respect to changes in microstructure. A governing 
differential equation and boundary conditions are obtained for the microstructural changes by posing 
the problem in a variational setting. The Eshelby stress tensor, a separate configurational stress, and 
thermodynamic driving (material) forces arise in this formulation, which is recognized as describing a 
process of self-assembly. An example is presented to illustrate the theoretical framework. 



1 Introduction 

The development of bioloKic al tissue consi sts of distinct processes of growth, remodelling and morphogenesis — 
a classification suggested bv lTaber In our treatment of the prob em, growth is defined as the addition 



or depletion of mass through processes of transport and reaction coupled with mechanics. As a result there 
is an evolution of the concentrations of the various species that make up the tissue. Nominally, these include 
the solid phase (cells and extra cellular matrix), the fluid phase (interstitial fluid), various amino acids, 
enzymes, nutrients, and byproducts of reactions between them. The stress and deformation state of the 
tissue also evolve due to mechanical loads, and the coupling between transport, reaction and mechanics. 

Remodelling is the process of microstructural reconfiguration within the tissue. It can be viewed as an 
evolution of the reference configuration to a "remodelled" configuration. While it usually occurs simultane- 
ously with growth, it is an independent process. For the purpose of conceptual clarity we will ignore growth 
in this paper, and focus upon a continuum mechanical treatment of remodelling. 

The microstructural reconfiguration that underlies remodelling is a motion of material point s in material 
space. An example of remodelling driven by stress is provided by the micrographs of Figuren]from 
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Figure 1: (a) Micrograph taken 3 days after plating of cells, shows a random distribution and orientation 
of tendon fibroblast cells in engineered tendon, (b) As growth occurs the cells organize into a more ordered 
microstructurc seen in this micrograph taken about a month after plating of cells. The horizontal alignment 
of cells corresponds to the orientation of a uniaxial stress that was imposed externally on the growing 
tendon construct. The alignment of cells along the stress axis is evidence of remodelling due to stress in the 
engineered tendon during growth. 



Remodelling also can be driven by the local density of the tissue's solid or fluid phases, availability of 
various chemical factors, temperature, etc. We assume that it is possible (through experiments or other 
approaches) to define a phenomenological law that specifies this evolution. Since any conditions that could 
drive remodelling vary pointwise through the tissue, such a "preferred" remodelled state will, in general, 
be globally incompatible. However, in its final remodelled state, the tissue is virtually always free of such 
incompatibilities (see Figure^. We therefore propose that a further, compatibilty-rcstoring material motion 
occurs, carrying material points to the remodelled configuration. 



Taber and HumnhrevI l|200l[) and lAmbrosi and Mollical 1)20021) have previously referred to remodelling 



However, the treatments in these papers are based upon concentration (or density) changes in growing tissue 
and the mechanics — mainly internal stress — that is associated with them . By our definitions, the se papers 
describe growth rather than remodelling. In a largely descriptive paper, iHumphrev and Raiaeopa l (2002J 
have proposed the evolution of "natural configurations" that seems closest to our ideas. To our knowledge, 
however, no quantitative treatment exists paralleling the ideas of microstructural reconfiguration, material 
motion, material/configurational forces and their relation to remodelling, as described in the present paper. 



2 Variational formulation: Material motion and material forces 

Figure |5] depicts the kinematics associated with remodelling. The preferred remodelled state is given by a 
tangent map of material motion K^:il,o x [0,?^] i— > GLP, where GL^ is the space of 3 x 3 matrices. The 
reference configuration is f^o C M^. Our first assumption is that is given. Futur e comm unications will 



address the derivation of such an evolution law from our experiments (see lCalve et al.l l|2003f) for preliminary 
results) . Since K' is generally incompatible (see Section ^ , a further tangent map of material motion, 
K'^-.flo X [0,T] M.^, acts to render the tissue of interest, B, compatible in its remodelled configuration, 
. The point-to-point map k: flo x [0,T] i— > carries material points from fig to il^ C M'^, the remodelled 
configuration. It is a material motion, and its compatible tangent map, K = dn/dX satisfies K — K'^K^ . 
The placement of material points in fl^ is X* = K[X,t). Further deformation, brought about by the 
displacement, u* , carries material points from to the spatial configuration Qf The deformation gradient 
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is i^* = 1 + du* /dX* . In this initial treatment we do not consider any further decompositions of F* . The 
overall motion of a point is t) = ^^i^i t) + '^*{^* 1 1) ° '^(^i t), and the corresponding tangent map is 

F = 1 + dip/dX. It admits the multiplicative decomposition F = F*K'^K^ . To reiterate upon the foregoing 
distinction between the material motion and deformation components of the kinematics, we emphasize that 
u* is a displacement, while k is a motion in material space. The corresponding tangent maps are F* (a 
classical deformation gradient) and K (a material motion gradient). 




Figure 2: The kinematics of remodelling 



2.1 A variational formulation 

We consider the following energy functional: 

n[M*,K] := j ^*{F\K'',X*)dV* 



at 

- J f* ■iu* + K)dV* -J t* ■{U*+ K)dA*, (1) 

n; an* 

where ip* = t/j* {F* , K'^ , X*) is the stored energy function. Observe that ip* is assumed to depend upon 
the compatibility- restoring material motion, K"^, in addition to the usual dependence on F* . Furthermore, 
material heterogeneity is allowed. The body force per unit volume in f2( is /*, and the surface traction per 
unit area on (?ri( is i* . Since the total motion of a material point is n + u* , the potential energy of the 
external loads is as seen in the second and third terms. 

Recall that the Euler-Lagrange equations obtained by imposing equilibrium with respect to u* (sta- 
tionarity of 11 with respect to variations in u*) represent the quasistatic balance of linear momentum in 

n;. 

'^U[uI,k] = where It* = u* + e^M* (2) 

e=0 



d£ 



Div*P* + /* = 0, in n*; P* N* = t* on dn*; P* (3) 

oF 
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Observe that the definition of P* resembles the constitutive relation for the first Piola-Kirchhoff stress if fl^ 
were the reference configuration. 

A final assumption is that the tissue also reaches local equilibrium with respect to k (stationarity of 11 
with respect to variations in k). The variational statement is: 



— 0, where = k + e5n, and if is fixed. (4) 

e=0 



The calculations are lengthy, but entirely standard, and yield the following Euler-Lagrange equations: 
- Div 



- Oinfi* (5) 



dX* 

N* = Oonar^*. (6) 



Observe that the Eshelby stress ^*1~F* P* makes its appearance. Hereafter, it will be denoted £. The 
term K'^ is a thermodynamic driving quantity giving the change in stored energy, ■0*, corresponding 
to a change in configuration, K'^ . It is stress- like in its physical dimensions and tensorial form, and we 
therefore refer to it as a configurational stress. Hereafter we win write S* = . 

Remark 1: A distinct class of variations can be considered than those in Specifically, consider 

— n[u*,/«e] = 0, where = /t + £(5/4, and c^, = X + K + £(5/« + u*, (7) 

de £=0 

which is distinct from Q in that variations on the material motion result in variations on the final placement 
as well. In this case too, the Euler-Lagrange equations © and JHl are arrived at. This is an important 
property: The system of equations governing the evolution of the microstructural configuration must be 
independent of the particular class of variations considered. 

2.2 Restrictions from the dissipation inequality 

The dissipation inequality written per unit volume in the reference configuration takes the familiar form, 

r: (FF-i) - ^ {dei[K]r ) > 0, (8) 

where t is the Kirchhoff stress defined in fit- Observe that det[X]^* is the stored energy per unit volume in 
r^o- Using the multiplicative decomposition F ~ F*K'^K^, and the defining relation for the configurational 
stress S* = K'^ , standard manipulations result in the following equivalent form: 

tF*'^ - det[ii:]|^) F* - det[K] {£ + S*) : ik" K''^^ 



'dF* 



-det[K]e: (k'k'k^^) - det[K]Hn^ ■ k > 0. (9) 

Adopting the constitutive relation, r = det[K]-^^F* (this is consistent with the observation that P* — 
^prr is related to the first Piola-Kirchhoff stress with fij as the reference configuration), results in the reduced 
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dissipation inequality 



-detK (£ + S*) : (k'K'^"^) 



-dei[K]S: (k^K'K-^) - det[K] ^.k>0, (10) 

which places restrictions on the evolution law for JC', and on the functional dependencies ?/;*(•, JC^, •) 
through S*, and 



3 Remodelling of one-dimensional bars 



In general, the examples of remodelling encountered in soft and hard biological tissue involve complex 
microstructural changes. Evolution laws for and the functional form, ij;* {• , K'^ , •) , to model these 
complexities are critical components for the successful application of the theoretical framework outlined in 
this paper. In future communications we will describe our experimental program to extract such constitutive 
information. However, the working of the formulation can be demonstrated by academic, but illuminating, 
examples. In the interest of brevity we restrict ourselves to a single example in this paper. 




L' / 



Figure 3: One-dimensional remodelling of bars. 



Consider two parallel bars, that may represent adjacent strips of a long bone (Figure |3J). We wish to 
consider a scenario in which each bar undergoes a preferred material motion to change its length from L to 
LI, i = 1,2. In general, this configuration is incompatible as the bars can attain different lengths. If they 
are required to remain of the same length in the remodelled configuration, further material motion occurs, 
resulting in a length L* for each bar. This is the remodelled configuration, Q^, in which the total material 
motion of each bar is k = L* — L. If the remodelling takes place under an external load, T, the bars each 
stretch to a final length I. The deformation is u* — I ~ L*. We examine the equations that govern the 
deformation and material motion by considering the following energy functional: 

n[u*, k] = ifc*(K + L - L{)^ + ^k*{K + L- U^f + 2 • ifcu*^ - r(M* + k), (11) 

where k* and k are spring constants for the material motion- and stretch-dependent portions of the stored 
energy, respectively. These portions are assumed to be separable. The theory of Section |21 results in the 
following relations: 

— = ^ 2ku* = T, (12) 
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an k ^ f L\ + Li\ 

In p2(l the standard relation is seen for the stretch of a hnear spring with effective stiffness 2k. The more 
interesting result is Observe that when L = ^{L\ + L2), material motion can occur, driven by stress, 

since u* = T/2k from H12() . In this case remodelling can be incompatible if ^ L\. However, remodelling 
does not drive material motion, k in this case. Instead there is stress-driven remodelling as described 
in Section ^ On the other hand, in the absence of an external load, material motion is obtained when 
L ^ \{L\ + U2). In this case the compatibility- restoring remodelling, motivated in Section ^^nd described 
by the tangent map, , also leads to overall material motion, k. 



4 Discussion and conclusion 

This paper has presented a theoretical framework for remodelling in biological tissue, where this phenomenon 
is understood as an evolution of the microstructural configuration of the material. The assumption that the 
material attains local equilibrium with respect to the evolution of its microstructure results in Euler-Lagrange 
relations in the form of a governing partial differential equation and boundary conditions. The final form of 
the equations and the results themselves depend critically upon the specified constitutive relations for the 
preferred remodelled state of the material, and the dependence of the stored energy upon the compatibility- 
restoring component of remodelling. These are open problems, and will be addressed in future papers by 
our group. The following points are noteworthy at this stage of the development of the theory: 

• This work does not deal with the approach to local equilibrium, which may take time on the order 
of days in biological tissue. Furthermore, the equilibrium state, being defined by the external loads, 
evolves upon perturbation of these conditions: Additional remodelling occurs in biological tissue when 
the load is altered. 

• The energy functionals in ^ and generalize to the Gibbs free energy of the body under constant 
loads and isothermal conditions. Further contributions that drive the process, such as chemistry or 
electrical stimuli, can be encompassed by the Gibbs energy. In such a setting the process we have 
described here would be termed a "self-assembly" in the realms of materials science or physics. 
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